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ABSTRACT 

We investigate the population of cool, evolved stars in the Local Group dwarf ellip- 
tical galaxy M32, using Infrared Array Camera observations from the Spitzer Space 
Telescope. We construct deep mid-infrared colour-magnitude diagrams for the resolved 
stellar populations within 3.5 arcmin of M32’s centre, and identify those stars that 
exhibit infrared excess. Our data is dominated by a population of luminous, dust- 
producing stars on the asymptotic giant branch (AGB) and extend to approximately 
3 mag below the AGB tip. We detect for the hrst time a sizeable population of ‘extreme’ 
AGB stars, highly enshrouded by circumstellar dust and likely completely obscured at 
optical wavelengths. The total dust-injection rate from the extreme AGB candidates is 
measured to be 7.5 x 10“^ Mq yr“^, corresponding to a gas mass-loss rate of 1.5 x 10“^ 
Mq yr“^. These extreme stars may be indicative of an extended star-formation epoch 
between 0.2 and 5 Gyr ago. 

Key words: galaxies: individual (NGC 205, M32) - stars: AGB - infrared: stars - 
stars: late-type - galaxies: stellar content 


1 INTRODUCTION 

Messier 32 (M32) is a dwarf elliptical galaxy in the Local 
Group and satellite of M31. Our relative proximity to M32 
(distance 785 ± 25 kpc; McConnachie et al. 2005) allows us 
to detect individual stars. Consequently, this provides us 
with the nearest opportunity to study the stellar popula- 
tions of an elliptical galaxy in detail. M32’s stellar pop- 
ulations have been extensively studied across many wave- 
lengths (e.g. Baade 1944; Rose 1985; Grillmair et al. 1996; 
Worthey et al. 2004; Monachesi et al. 2011; Sarajedini et al. 
2012). However, these populations have not been charac- 
terised in the mid- infrared (mid-IR). 

It is thought that M32 has two main populations: an 
intermediate-age, metal-rich ([Fe/H] = ~-|-0.1) population 
thought to be ~2-8 Gyr old (Rose et al. 2005; Coelho et al. 
2009), and an old (8-10 Gyr) stellar population with slightly 
sub-solar metallicity ([Fe/H] = -0.2; Monachesi et al. 2011). 
Additionally, a small number of ancient, metal-poor ([Fe/H] 
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= -1.42) RR Lyrae variables in M32 have recently been 
identified (Fiorentino et al. 2010, 2012), and a young pop- 
ulation of 0.5-2 Gyr stars is probably present at all radii 
(Monachesi et al. 2011). The extremely high surface bright- 
ness and sharply peaked brightness profile towards the nu- 
cleus of M32 indicates a strong increase in the stellar density. 
Given the extreme crowding, the stellar populations within 
the core (r < 1' of the nucleus) of M32 are incompletely 
probed. This high surface brightness also prevents studies 
of the fainter stars in M32. 

The unique nature of M32 in the Local Group and 
the controversies over its origins make M32 a tantalising 
prospect for stellar population studies (proposed forma- 
tion scenarios include a true elliptical galaxy or a severely 
truncated early-type spiral; Faber 1973; Bekki et al. 2001; 
Ghoi et al. 2002; Chilingarian et al. 2009; Kormendy et al. 
2009). Our immediate goal is to assess the infrared stellar 
populations of M32 using Spitzer /I'R.A.C observations and to 
study the circumstellar dust emission from late-stage AGB 
stars. 

The Infrared Array Camera (IRAC; Fazio et al. 2004b) 
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on-board the Spitzer Space Telescope (Werner et al. 2004) 
provides a unique window to probe M32’s cool, luminous 
stars, whose spectral energy distributions peak in the near- 
IR and which often have an infrared excess due to thermal 
emission from circumstellar dust. Indeed, the dustiest, ex- 
treme AGB stars are often obscured at optical and near- 
IR wavelengths, and thus can only be detected by their 
infrared emission (Meixner et al. 2006; Blum et al. 2006; 
Mould et al. 2008; Boyer et al. 2011; Gerke & Kochanek 
2013). 

‘Extreme’ AGB (x-AGB) stars experience a consider- 
able mass loss and are thought to dominate the mass return 
to the interstellar medium (ISM) in Local Group galaxies 
(Riebel et al. 2012; Boyer et al. 2012). These stars were first 
identified by Blum et al. (2006) as AGB stars whose near- 
IR stellar flux is heavily extinguishes by circumstellar dust, 
thus their K-band magnitude is often fainter than the tip of 
the RGB (see Boyer et al. 2011). The criterion Blum et al. 
(2006) used to select x-AGB star candidates was bases on 
their location in the [8.0] versus J — [8.0] colour magnitude 
diagram, where a break is seen in the carbon-star sequence 
at J — [8.0] = 3.5. Later this selection criteria was modified 
by Boyer et al. (2011) to focus on stars brighter than the 
tip of the RGB at 3.6 /un with [3.6]-[8.0]> 0.8; enabling 
the identification of objects completely obscured at optical 
and near-IR wavelength. In this paper we adopt the latter 
selection criteria for x-AGB stars (see Section 4.3). 

By definition, the x-AGB class represents the high mass- 
loss rate evolved stellar population in a system. Identifying 
these sources is important for constraining the dust budget 
of a galaxy (Riebel et al. 2012). Particularly as the most- 
extreme sources: AGB stars experiencing the final ‘super- 
wind’, only constitute a small faction of the x-AGB class. 

In this paper we describe the results of Spitzer IRAG 
and MIPS observations of M32. Our observational data, re- 
duction and photometry is described in Section 2, and we 
discuss contamination by M31 and foreground stars in Sec- 
tion 3. Infrared colour-magnitude diagrams (CMDs) and lu- 
minosity functions are presented in Section 4; we use these 
to determine the integrated dust return to the ISM (Sec- 
tion 4.4). Finally the conclusions are summarised in Section 
5. 


2 OBSERVATIONAL DATA 
2.1 Observations 

Spitzer observations of M32 (Program Identification [PID] 
3400, PL M. Rich) were obtained on 18 January 2005 UT 
with IRAG (Fazio et al. 2004b) at 3.6, 4.5, 5.8 and 8 /rm. 
For each IRAG channel, 92 (5.2’ x 5.2’) frames at 23 dither 
positions were taken in a cycling pattern with 30-second 
exposures to build redundancy against outliers and arti- 
facts, resulting in total integration time on M32 of 2760s 
over most of the map. Observations were centred at R.A. 
= 00'*42’"41L6, deck = -f40'’55'53".6 (J2000.0) and cover 
an area of approximately 6'x7' in all four IRAG channels 
around the centre of M32, plus an off-held region (offset by 
one pointing) of the same size to the north-west for IRAG 
bands 2 and 4 or south-east for IRAG bands 1 and 3. Fig- 
ure 1 shows this field of view overlaid on the MIPS 24-^m 



Figure 1. Location of our Spitzer IRAC pointings towards M32 
superimposed on a MIPS 24 pra mosaic (Gordon et al. 2006; 
Bendo et al. 2012). The region observed with IRAC [3.6] and [5.8] 
is shown in red and the IRAC [4.5] and [8] observations are shown 
in blue. 


mosaic of M31. No dedicated background control fields were 
observed to complement these observations. 


2.2 Data Reduction 

The raw data was processed by the Spitzer Science Center 
(SSC) reduction pipeline version S18.7.0; this removes elec- 
tronic bias, subtracts a dark image, applies a flat-field cor- 
rection, and linearises the pixel response. The resulting Basic 
Calibrated Data (BCD) images were corrected for various in- 
strumental artifacts (e.g. muxbleeds and column pull-down 
effects) and combined using the mopex reduction package 
(Makovoz & Marleau 2005) to produce a single mosaic for 
each channel. 

The MOPEX OVERLAP routine was implemented to 
match the backgrounds of individual frames in overlap- 
ping areas of the images producing a smooth background. 
The MOSAICKER, in addition to image interpolation and co- 
addition, eliminated cosmic rays and other outliers in the 
data. The final IRAC mosaics are not sub-sampled; thus 
each image has a pixel scale of 1.22" pixeP^ and a Point 
Source Function (PSF) Full With Half Maximum (FWHM) 
of ~1.9 arcsec. 

24 pm data from the Multiband Imaging Photome- 
ter for Spitzer (MIPS) were used in part of this analysis. 
The data, which was originally published by Gordon et al. 
(2006), were reprocessed using the MIPS Data Analysis 
Tools (Gordon et al. 2005) and additional data processing 
steps described by Bendo et al. (2012). The image has a 
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Table 1. Catalogue of Point Sources for M32. 


Apparent Magnitude 


SSTM32 # 

Source ID 

RA 

Dec 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

0001 

J004235.04-I-404839.2 

10.6460 

40.8109 

13.885 ± 0.051 

12.694 ± 0.088 

12.243 ± 0.065 

11.446 ± 0.063 

0002 

J004227.59-I-405300.2 

10.6150 

40.8834 

15.378 ± 0.241 

15.043 ± 0.215 

999.99 ± 999.9 

15.317 ± 0.125 

0003 

J004238.90-I-404830.5 

10.6621 

40.8085 

15.004 ± 0.184 

13.739 ± 0.116 

13.347 ± 0.043 

12.543 ± 0.055 

0004 

J004241. 59-1-404917.7 

10.6733 

40.8216 

17.350 ± 0.242 

16.754 ± 0.208 

999.99 ± 999.9 

999.99 ± 999.9 


Notes: The source ID follows the standard Spitzer naming convention, giving the truncated (J2000) coordinates. 
Only the first entries are shown for guidance, a full version is available in the electronic edition. 



Figure 2. Three-colour image of M32. Blue is 3.6 /rm, green is 
8.0 /rm and red is 24 /rm. Point sources with a dust excess become 
visible at 8.0 /rm and 24 fira and appear red in the image. 


pixel size of 1.5 arcsec, a PSF FWHM of 6 arcsec, and a cal- 
ibration nncertainty of 4 per cent (Engelbracht et al. 2007). 

Initial inspection of the images revealed a population of 
clearly resolved point sources at 8 /tin, associated with M32 
but not obviously matched at shorter wavelengths. Figure 2 
presents a three-colour image combining the 3.6, 8.0, and 24 
/im observations; the red point sources are clearly visible. 
Additionally, we see no emission indicative of interstellar 
dust within M32. This would be particularly evident in the 
8.0-/tm map, which is sensitive to non-stellar emission from 
polycyclic aromatic hydrocarbons (PAHs), or the MIPS 24- 
/im mosaic, which is sensitive to cool dust. 


2.3 Photometry 

Point sources were extracted from individual IRAC frames 
and on the mosaicked images using PSF fitting with the 
DAOphot II and allstar photometry packages (Stetson 
1987). To determine the shape of the PSF, DAOphot re- 
quires selection of isolated stars in the field; this is quite 
challenging in very crowded fields such as M32. PSFs were 


Table 2. Number of point sources detected in M32 


Number 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

Total 

All detections 

2177 

2457 

380 

575 

5589 

Final Catalogue 

1270 

1408 

349 

468 

1552 


created from a minimum of 10 isolated bright stars in each 
IRAC channel, and sources 4cr above the local background 
were chosen for extraction. 

Allstar was then used to fit PSFs to all sources found 
on the image. This is an iterative process that simultane- 
ously fits each source in the frame with the PSF profile 
and subtracts converged sources from the input image. As 
sources are removed, new estimates for the local background 
are calculated, improving the flux estimates for crowded and 
faint sources. Sharpness and roundness cuts were employed 
to eliminate from the sample extended sources, cosmic rays 
and unrecognised blends that are broader or narrower than 
the PSF. Cosmic rays are also eliminated from the point- 
source list when detections in the same band are combined. 
The resulting source lists are then cross-correlated with the 
point-source photometry from the mosaicked data. 

The flux densities and uncertainties are colour- 
corrected using a 5000 K blackbody, according to the 
method described in the IRAC Data Handbook, ver- 
sion 3.0^. Additionally, a pixel-phase-dependent correction 
(Reach et al. 2005) was applied to the 3.6-/tm photome- 
try. Finally, magnitudes relative to Vega were derived using 
the zero- magnitude flux calibrations provided in the Spitzer 
IRAC Data Handbook. Figure 3 shows the representative 
photometric uncertainty as a function of source magnitude. 
The IRAC instrument also has an absolute calibration un- 
certainty of 3 per cent (Reach et al. 2005), although we do 
not include this uncertainty in our photometric errors. 

Table 1 lists the 1552 sources included in the final point- 
source catalogue. For a source to be included in the cata- 
logue we require it to be detected in at least two IRAC chan- 
nels; we consider a source to be a match if their centroids 
are within a T'radius. The point sources detection statistics 
are summarised in Table 2. 

Completeness limits for each wavelength were deter- 
mined via the extraction of artificial inserted stars. These 
false stars were placed at random pixel locations and have a 
limiting magnitude ~ 2 mag fainter than the faintest objects 


1 


http: / /ssc . spitzer . caltech.edu/irac/dh/iracdatahandbook3 . 0 .pdf 
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Figure 3. Photometric uncertainty and completeness fraction as a function of apparent magnitude for each IRAC band. The one- 
dimensional completeness fraction as a function of magnitude is shown in the bottom panels. Limits are averaged over the galaxy: 
photometry is more complete in the outskirts of the galaxy compared to the core. 


recovered from our images. To avoid a significant increase 
in the source density of the field we limit the fraction of 
artificial stars introduced per image to ~ 5 per cent of the 
total point sources detected. PSF-fitting photometry was 
then performed on the modified image. This process was re- 
peated multiple times for each image to statistically analyse 
the photometric completeness of each band. Artificial stars 
are considered to be recovered if they are within a one-pixel 
radius of the input position and their magnitude differs by 
\5m\ ^ 1 from their input magnitudes. If the magnitude dif- 
ference is greater than this we consider the sources to be 
an unresolved blend of two or more stars and hence we do 
not consider them as recovered. The bottom panel in Fig- 
ure 3 compares the fraction of injected and recovered stars 
for each IRAC band. 

The different stellar density conditions found across our 
field of view have a large impact on the fraction of point 
sources recovered. Severe stellar crowding in the central re- 


gions of M32 results in a high degree of blending; this proves 
to be the major source of incompleteness. Furthermore, for 
[3.6] and [4.5] the completeness factor of recovered stars in a 
given magnitude bin is much lower than at [5.8] and [8.0]. At 
shorter wavelengths, the increased sensitivity and brighter 
stellar emission results in a higher detection rate and thus a 
higher fraction of blended sources in crowded regions. 


In regions of high surface brightness, unresolved blended 
sources in our catalogue will show a magnitude enhancement 
compared to individual point sources. Figure 4 plots the ap- 
parent magnitude as a function of radial distance. The rel- 
atively flat profile indicates that blended sources form only 
a minor component of our catalogue. Any sources in the in- 
ner regions that show a significant magnitude enhancement 
(> 5 (t above the mean at R > 1.5 ') are removed from our 
catalogue. 
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Figure 4. Distribution of the apparent magnitude with radial 
distance for all point sources detected at 8.0 /rm. Towards the 
nucleus of M32 crowding is severe; sources with a magnitude en- 
hancements in the core regions are probable blends. 
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3 CORRECTING STAR COUNT NUMBERS 
3.1 M32 field contamination 

Typically, Galactic foreground stars as well as unresolved 
background galaxies will contribute to the source contam- 
ination in our catalogue. These contaminants will have 
a spatially uniform distribution across our field of view. 
To estimate the degree of contamination from foreground 
sources we use the trilegal population synthesis code of 
Girardi et al. (2005) and we extrapolate the extragalactic 
background contamination from Fazio et al. (2004a). In our 
IRAC field, we find that the contamination from foreground 
Galactic stars accounts for ~ 6 per cent of the point sources 
detected at 3.6 /rm in our photometric catalogue. Compara- 
tively, the upper limit for the number of background galaxies 
is ~ 18 per cent of objects detected at 8 /rm. However, the fi- 
nal number of background galaxies included in our catalogue 
may be considerably lower as our point source extraction 
method rejects extended non-stellar objects. An alternative 
estimate of the foreground star and background galaxy con- 
tamination may be obtained via colour-magnitude diagrams 
(Section 4.3). 

M32 has a projected separation of 24' (5.4 kpc) from 
M31’s centre and is seen against the disc of M31. Disc and 
halo stars from M31 are the dominant source of contam- 
ination in our M32 field. Furthermore, stars belonging to 
the nearby M31 galaxy will have similar IR colours to the 
evolved stellar population of M32. This makes it extremely 
difficult to isolate the individual stars belonging to M32, 
and the presence of some intruders is unavoidable. Com- 
bined with crowding, contamination from M31 is the most 
important limitation in the analysis of M32’s rich stellar 
field. 

As contamination from M31 will not be homogeneously 
distributed across the field of view, the best estimate of the 
fraction of contaminating objects would be from a suitable 


Figure 5. Source density profiles for M32. Spitzer resolved the 
bright stellar populations in the outer halo of M32; however, the 
extremely high surface brightness and crowding in the bulge hin- 
ders detections in the central regions. The dashed line marks the 
radial limit where we have coverage in all four IRAC bands. 


control field located at the same isophotal level in the outer 
disc of M31 but away from M32’s nucleus. Unfortunately, 
archival Spitzer IRAC observations of M31 and its surround- 
ings only mapped the outer regions of M31 for about 120 
seconds per sky position (Barmby et al. 2006; Mould et al. 
2008). These observations do not provide sufficient depth 
to serve as a control field. In the following subsections, we 
describe the methods taken to mitigate the impact of M31 
disc stars on our sample. 


3.2 Source density profiles 

Using the radial density profiles we attempt to measure the 
IR contamination from M31 and Galactic foreground stars 
in each field of view. To measure the stellar radial density 
profile of M32 we counted the sources contained in concen- 
tric annuli from its centre (where the bulk stellar population 
peaks) in steps of 0.5'. The number of sources per unit area 
was then obtained by dividing this value by the correspond- 
ing area covered by the annuli after correcting for regions 
not covered by our observations. Figure 5 shows the source 
density profiles for M32. 

In most channels the source density declines signifi- 
cantly toward the centre of M32, where incompleteness due 
to crowding effects in the core of M32 becomes significant. 
Both the 5.8- and 8-fim bands show a smooth decline in 
source density beyond a radius of R > 1.0', until contam- 
ination from field stars and M31 begins to dominate the 
source density. This is particularly noticeable in the 8-^m 
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Figure 6. King profiles fit to the background corrected IRAC 5.8 
and 8 micron density profiles for M32. The data deviates from the 
King profiles due to source confusion inside 1 arc-minute. This 
indicates that we are underestimating the stellar density in the 
centre of M32. 


band beyond R > 4.5', where there is a steady rise in source 
counts towards M31. 

The flattening of the star counts at 7? = 2.5 - 4.0' was 
used to estimate the contamination from M31 and fore- 
ground stars. This corresponds to a point-source density of 
~5 sources per arcmin^ for the 5.8 /rm band and -^7 sources 
per arcmin^ for the 8 /rm band. The flat stellar density pro- 
files in the IRAC 3.6 and 4.5 fim bands beyond R > 2.0' in- 
dicates a high degree of point source confusion with M31. 
Because of these high stellar surface densities we exclude 
these channels from the M32 stellar population analysis in 
Section 4.1. 

To determine the radius at which our data becomes con- 
fusion limited we fitted the radial profiles of M32 with an 
empirical King (1962) profile of the form: 


Table 3. King profile function best-fit parameters 


Parameter IRAC 5.8 fim IRAC 8.0 pm 


lo 

Rc (arcmin) 
Rt (arcmin) 


1103.68 

0.241 

3.36 


1236.96 

0.245 

3.16 


2.5'), uncertainties in the background source density causes 
deviations from the profile shape. This also limits how well 
we can constrain the tidal radii. Table 3 lists the best-fit 
parameters. 

Due to the flat stellar-density profiles in the IRAC 3.6- 
and 4.5-/rm bands we were unable to determine the corre- 
sponding King profiles for these channels. These flat profiles 
arise because of the low number of reliable detections within 
90 arcseconds of M32’s core and contamination from M31’s 
disc at large radii where there is limited radial coverage for 
all angles. 

Based on the integrated surface brightness profile at 
optical wavelengths, Choi et al. (2002) find that the inner 
regions of M32 (r < 150"; 0.6 kpc) are predominantly spher- 
ical and well characterised by a de Vaucouleurs profile, while 
the outer regions show evidence of a residual elongated disc. 
We find a similar break at r « 150", however, at large radii 
the number density of point sources per unit area becomes 
progressively more uncertain, and so our star-count analysis 
is of limited use for assessing a potential disc population. 

At M32’s distance, the tidal radii in both the 5.8- and 
8-pm bands correspond to a projected distance of 0.7 kpc. 
This value is somewhat lower than that given by Choi et al. 
(2002) of Rt = 310" (1.2 kpc), suggesting that the mid- 
IR stellar population maybe slightly more centrally concen- 
trated than the underlying unresolved optical sources. This 
would be consistent with a more centrally concentrated pop- 
ulation forming at the centre of M32’s gravitational poten- 
tial well, where gas would naturally collect. Such a distri- 
bution in the stellar population is consistent with the idea 
that the progenitor of M32 was a larger disc galaxy striped 
by interaction with M31. This scenario was proposed by 
Bekki et al. (2001) to explain M32’s relatively high metal- 
licity, its low gas content and the lack of a globular clus- 
ter system. However, it is also consistent with the bright 
mid-IR population arising in a more metal rich or younger 
population, both of which can be plausibly more centrally 
concentrated, without invoking mergers. 


HR) = lo 


(1 + 


)-o-5_(l + 


R < 


where I (R) is the density profile at projected radius R, Iq is 
a normalising factor, Rc the core radius, and Rt is the pho- 
tometric tidal cut-off radius. The best-fit King profiles to the 
background-corrected IRAC 5.8- and 8-fim density profiles 
are shown in Figure 6. In both channels the data deviate 
from a King profile and become confusion limited at i? < 1', 
which prevents us from studying the stellar populations in 
the innermost part of M32. At larger radial distances {R > 


3.3 Stellar spatial distribution 

Due to the limited spatial coverage, star counts measured 
in the off-field regions were also used to estimate the M32 
^*(^keld contamination. These regions (shown in Figure 1) were 
obtained in conjunction with our primary observation of 
M32 and thus have the same depth as the M32 field. Star 
counts were measured using concentric annuli to sample a 
sub-section of each field, and cover the same solid angle as 
those used to measure the source density profile of M32. The 
off-field annuli centred near a =0:42:56.10, 5 = -1-40:45:48.38 
sampled the outer disc of M31 where we have coverage at 
3.6 and 5.8 fim. The other set of off-field annuli was centred 
at Q =0:42:27.75, <5 = -1-40:57:48.76 in the field closer to the 
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main body of M31 observed at 4.5 and 8 /rm. In these re- 
gions the stellar density of M31 should average out across 
the field of view, producing a flat density profile. 

Located furthest away from M31, the IRAC 5.8-fim 
field suffered the least contamination, with an average back- 
ground point-source density of 2.5 sources per arcmin^, while 
the 8-^m field (which intersects the inner disc of M31) has 
an average density of 11 sources per arcmin^. Based on the 
ratios of sources in the off-fields, we expect < 20 per cent 
of the point sources in the on-field region to be contami- 
nants from M31 in the 5.8- and 8-/rm bands. This assumes 
that objects associated with M32 do not contribute to the 
number density in the off-held regions. 

4 RESULTS & DISCUSSION 
4.1 M32 luminosity functions 

In this section, we examine the mid-IR stellar luminosity 
functions of M32. These are useful in constraining the star- 
formation history and provide an estimate of the dust emis- 
sion. The luminosity functions for each of the IRAC bands 
were constructed by sorting the point sources into bins based 
on their magnitude. The bin size of the luminosity functions 
is set to 0.2 mag, to account for photometric errors in our 
catalogue. To reduce the effects of contamination by M31 as 
much as possible, while still retaining a reasonable sample 
size, only sources which are detected in at least two wave- 
bands and that are within 3.5' of M32’s centre were counted. 
It should also be noted that the luminosity function at 3.6 
^m will be suppressed due to the extreme crowding in the 
centre of the galaxy. Furthermore, if the centre contains a 
different stellar population than the outer part of the galaxy, 
the shape may also be affected. 

Comparisons of the 3.6-^m luminosity function with 
and without stars in the central 1' yield similar results; this 
is not unexpected as the increase in crowding near the core 
limits the number of stars that can be reliably resolved. At 8 
^m the situation is different. At these wavelengths crowding 
has less of an impact and we are sensitive to sources with 
r > 0.6' of the core. Removing sources within r < 1' effec- 
tively scales the luminosity function downwards and we see 
a very slight change in the gradient. At the faint end of the 
luminosity function we detect very few sources within 1'; it 
is here where the effects of crowding and the high surface 
brightness has some impact. 

The uncertainties in the luminosity functions are set 
both by Poisson statistics of the sample and by the un- 
certainties in the background counts for each magnitude 
bin. Defining a complete, uncontaminated sample consisting 
solely of M32 sources is not possible with currently available 
data; thus for the present survey, the luminosity functions 
presented here provide an upper limit to the true luminosity 
distribution, as some contamination may be present in each 
magnitude bin. 

Figure 7 shows the Spitzer IRAC luminosity functions 
at 3.6 and 8.0 fim for M32. The profile shape at 3.6 pm is 
similar to the variable stars in M33 (McQuinn et al. 2007), 
but it has a steeper falloff at bright magnitudes compared 
to all detections in the Magellanic Clouds. This is likely due 
to a difference in the metallicity of the stellar populations 


and hence the ratio between C-rich and M-type AGB stars 
(Mouhcine & Langon 2003; Battinelli & Demers 2005; Cioni 
2009). 

At 8.0 pm the luminosity function drops smoothly with 
magnitude and is similar in appearance to both the Mag- 
ellanic Clouds and M33; at this wavelength the luminos- 
ity function probes the dust emission rather than the pho- 
tospheric temperature. The large number of point sources 
in M32 indicates a significant population of dust-producing 
stars. 

The tip of the AGB can be estimated from the lumi- 
nosity functions by identifying the magnitude where the 
source count decreases significantly. Although source counts 
are low, an inspection of the luminosity functions suggests 
that the tip of the AGB for M32 is located near Ma.g 
— 10.8± 0.6 mag. M32 candidate stars brighter than the 
tip of the AGB are either luminous RSGs in M31 or fore- 
ground sources. The position of the tip of the AGB indicates 
that the IR-stellar population of M32 appears older than 
the Magellanic Clouds, which have experienced an extended 
star forming history. The location of the tip of the red giant 
branch (TRGB) is relatively insensitive to the metallicity 
of the host population (Cioni & Habing 2003), and lies at 
AI3.6 ~ —6 mag, below the completeness limit of our survey. 
We therefore assume that the majority of the point sources 
detected in our catalogue are intermediate-age AGB stars, 
with some small contamination by other stellar types. 

In contrast to the luminosity function of globular clus- 
ters where there is only a sparse population of AGB stars 
above the TRGB (Figure 8), M32’s luminosity functions are 
strongly dominated by a bright evolved stellar population. 
Stars in Galactic globular clusters are typically found to be 
10-14 Gyr old (e.g. Roediger et al. 2014), while bright car- 
bon stars are predicted to be much younger; with ages in 
the region of 8 Gyr to 100 Myr, depending on the mass 
and metallicity of the system (Karakas & Lattanzio 2007; 
Marigo & Girardi 2007). The majority of stars detected at 
8 pm in M32 are thus much younger than the Population II 
stars in globular clusters, and the numbers suggest we are 
tracing a metal-rich population of intermediate age. 

We detect a large number of dusty stars compared to 
the Local Group dwarf irregular galaxies WLM, IC 1613, 
Phoenix, LGS 3, DDO 210, Leo A, Pegasus dirr, and Sex- 
tans A (Jackson et al. 2007a, b; Boyer et al. 2009) commen- 
surate with the greater luminosity of M32. These composite 
stellar populations span a wide range in age and metallicity, 
however, any further comparison is limited due to the low 
numbers of stars in these systems. 

4.2 IRAC Colour distribution 

Figures 9 and 10 show histograms of colours for point sources 
in our catalogue. Oxygen-rich giants typically have a slightly 
blue [3.6]-[4.5] colour distribution, due to photospheric SiO 
and CO absorption at 4.08 pm and 4.66 pm, respectively, 
while HCN and C2H2 absorption features at 3.5 pm and 
3.8 pm, within the 3.6 pm band cause carbon stars to ap- 
pear red (Marengo et al. 2007). Foreground main-sequence 
dwarfs are also slightly reddened. In Figure 9 the blue- 
shaded area shows all the sources detected at 3.6 and 4.5 pm; 
this population is a mix of foreground dwarf stars, contami- 
nants from M31 and AGB stars in M32. Evolved stars under- 
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Absolute Magnitude 


Figure 7. Comparison of the M32, M33 and Magellanic Cloud luminosity functions. The dotted-dashed line marks the 50% completeness 
limit. The dotted line marks the 3.6 /rm TRGB for the LMC and the dashed line shows the TRGB for the SMC. 


going significant mass loss become dust-enshrouded and ex- 
hibit a strong mid-IR excess. Consequently, objects detected 
at [8.0] can be tentatively separated into oxygen-rich and 
carbon stars; this restriction helps to minimise contamina- 
tion from foreground sources. As the IRAC sensitivity limit 
for the 8-pm data is brighter than at 3.6 and 4.5 /rm: bright 
evolved AGB stars with strong circumstellar dust emission 
are preferentially detected over fainter, less-evolved AGB 
stars. This favours carbon stars and might explain the low 
number of oxygen-rich objects in our sample. 

Figure 10 separates the [3.6]-[8.0] colour distribution as 
a function of magnitude, following McQuinn et al. (2007). 
The bright sources with M 3.6 < —9.5 have a binomial distri- 
bution with peaks at [3.6] — [8.0] « 0 (which corresponds to 
AGB stars without dust) and [3.6] — [8.0] « 1.5 (sources red- 
dened by circumstellar dust emission). Based on the colour 
distribution of variable stars in M33, carbon star candidates 
are thought to fall in the range 0.5 < [3.6] — [8.0] < 2.5. 
Other non-variable sources reddened by circumstellar dust 


emission, for example, young stellar objects would have 
[3.6] - [8.0] > 2.8. 

4.3 IRAC Colour-Magnitude Diagrams 

Mid-IR colour-magnitude diagrams (GMDs) for all cata- 
logue point sources are presented in Figures 11 and 12. At 
these wavelengths stellar temperature no longer affects the 
GMDs (as the SED is described by the Rayleigh- Jeans tail), 
and any features seen are due to the absorption and emis- 
sion by circumstellar molecular and dust species. The well- 
defined vertical branch with zero colour (]3.6]-]8.0] « 0) in 
these diagrams traces foreground stars and dust-free blue- 
supergiants in M31. The relatively low number of sources 
that fall on this branch indicates a nominal level of contam- 
ination from foreground sources within our field of view. 

Based on the maximum luminosity at 3.6 /rm and the 
slightly red vertical [3.6]-]4.5] colour distribution, the fore- 
ground sources are probably main-sequence dwarfs. Unlike 
giants, dwarf stars have water vapour absorption features 
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Figure 8. Comparison of the luminosity functions of M32 and 47 Tuc. The red line represents all sources detected within 3.5'of M32’s 
centre; the black line represents the IRAC deep ‘Rood’ observations of 47 Tuc (Origlia et al. 2007; Boyer et al. 2010) and the green line 
represents the SAGE-SMC observations of 47 Tuc (Boyer et al. 2010). This indicates that the recovered stars in M32 are significantly 
brighter than the TRGB. 


at 3.6 /tm and lack CO absorption at 4.5 /tm resulting in 
a shift towards the red (Merrill & Stein 1976); foreground 
giant stars would also appear brighter. 

As we probe to fainter magnitudes, photometric errors 
become more significant. This may be responsible for the 
large spread in colour for [3.6] ^ 16.5 in the [3.6] vs. [3.6]- 
[4.5] colour-magnitude diagram (Figure 11, top panel). It 
should be noted that sources towards the centre of M32 
may be not accurately extracted at one or more wavelengths 
due to crowding. Although care had been taken to remove 
blended sources from the catalogue, blends may potentially 
affect the derived colours by enhancing the IRAC flux. We 
expect this to influence < 4 per cent of the sources in the 
CMDs, if we assume all detections with magnitude enhance- 
ments greater than 3cr are blends. 

The [8.0] vs. [5.8]-[8.0] colour-magnitude diagram (Fig- 
ure 12, bottom panel) provides the best representation of 


the ‘true’ stellar population of M32, as these bands are least 
affected by contamination from M31. 

The stellar sequences traced by the diagonal branches 
in the CMDs are difficult to isolate, making it harder to ac- 
curately classify the different populations of cool evolved 
stars. Using the mid-IR colour-classification scheme from 
Boyer et al. (2011), we identify a large population of dust- 
enshrouded evolved stars (x-AGB stars) with [3.6] — [8] > 
0.8. This corresponds to 55 per cent of the sources with 
detections in both channels. In non-star-forming galaxies 
like M32 contamination from young stellar objects should 
be negligible; however unresolved background galaxies may 
contaminate the x-AGB sample. Contamination from un- 
resolved background galaxies may be reduced by requiring 
x-AGB sources to be brighter than an empirical boundary 
set for extra-galactic objects. Thus we apply an additional 
colour-cut in the [3.6] — [8] vs. [8] GMD, which is adjusted for 
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Figure 9. The IR.AC [3.6] — [4.5] colour distribution of all detected 
sources in M32, shown in blue. Objects which have a measured 
[8.0] mag in M32 are denoted by red. Stars detected at [8.0] and 
appearing to the right of the dashed line are probably carbon-rich. 



[3.6]-[8] 

Figure 10. The IRAC [3.6] — [8.0] colour distribution of sources 
in M32. Bright sources with M3.6 < —9.5 are shown in blue while 
faint sources (M3.6 > —9.5) are shown in red. The dashes vertical 
lines mark the range used to select Carbon star candidates. 


M32’s distance from the classification by Boyer et al. (2011) 
and defined as: 

[8] = 17.4 - 1.33 X ([3.6] - [8]); [3.6] - [8] < 3.0. (2) 

There are 11 sources that do not fall within this range, re- 
ducing the final number of x-AGB candidates to 110. The 
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Figure 11. Spitzer IRAC CMDs of M32. For all panels the y-axis 
is the apparent 3.6 pm magnitude. The red triangles indicate a 
8-pm counterpart. The initial colour-cut used to select x-AGB 
candidates is indicated by the vertical dashed line, while the hor- 
izontal line represents the theoretical TRGB for M32. The error 
bars show representative uncertainties; they do not vary signifi- 
cantly with colour. 


colour-cuts used to identify the x-AGB candidates are plot- 
ted in Figures 11 and 12. Although an effort was made to 
remove contaminating sources from the final list of the x- 
AGB candidates, we still expect a moderate contamination 
from background galaxies. Other colour-cuts used to iden- 
tify 0-AGB and G-AGB stars require additional near-IR 
photometry. 

The x-AGB sample is thought to be dominated by 
carbon stars (van Loon 2006; Boyer et al. 2012), although 
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Figure 12. Spitzer IRAC CMDs of M32. For all panels the y-axis 
is the apparent 8.0 pm magnitude. The colour-cuts used to select 
x-AGB candidates are indicated by the dash-dot lines lines. The 
error bars show representative uncertainties; they do not vary 
significantly with colour. 
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some cross-contamination by oxygen-rich AGB stars is ex- 
pected. At LMC and SMC metallicities, 97 per cent of the 
x-AGB candidates were matched to a carbon-rich GRAMS 
model (Riebel et al. 2012). Spectroscopic observations of the 
Magellanic Clouds confirm this: less that 10 per cent of 
the candidate x-AGB stars were identified as oxygen-rich 
stars (Woods et al. 2012; Jones 2013). Previous estimates of 
M32’s metallicity suggest that the galaxy is predominantly 
metal rich, with a near-solar [Fe/H] value which trails off to 
[Fe/H]~-1. Consequently, we expect to see a larger (> 4 per 
cent) fractional contribution from oxygen-rich AGB stars 


Figure 13. IRAC [3.6] — [4.5] vs. [5.8] — [8.0] colour-colour di- 
agram, showing the 139 point sources identified at all bands in 
our M32 catalogue. The extreme AGB stars in our sample are 
indicated by red stars. 


in the x-AGB class (Boyer et al. 2013), although without 
spectroscopic confirmation how much more is unclear. Due 
to sensitivity limits these 0-rich stars are likely to be more 
massive than 4 Mq and have undergone hot-bottom burning 
for the star to remain 0-rich. 

At metallicities comparable to M32, AGB stars will be- 
come carbon-rich if their initial mass is in the region of 
1.2 - 3.4 Mq (Straniero et al. 1997; Marigo & Girardi 2007; 
McDonald et al. 2012, figure 9); this corresponds to a stellar 
age between 0.2 to 5 Gyr. This boundary is dependent on 
the choice of evolutionary model and may potentially extend 
to slightly higher ages (up to 8 Gyr). Thus the brightest x- 
AGB stars may be the progeny of a younger subset of the 
stellar population of M32. This agrees well with suggestions 
that the metal-rich component of M32 was formed during a 
period of star formation that occurred less than 7 Gyr ago 
(Monachesi et al. 2011). 

Since our data is more sensitive to carbon stars which 
trace a younger stellar population than the majority of 
oxygen- rich AGB stars (which have a lower initial mass), 
we cannot exclude the possibility that an older oxygen- 
rich population is also present in M32. Evidence for mul- 
tiple star-formation events in M32 has been presented by 
Monachesi et al. (2012) using stars near the main sequence 
turnoff. These multiple star-formation events may have been 
induced by tidal interactions with M31 or through accretion 
of gaseous material (Bekki et al. 2001; Block et al. 2006; 
Hammer et al. 2010). An extended period of star-forming 
activity associated with a much larger disc-dominated 
galaxy would increase the numbers of old AGB stars; unfor- 
tunately these faint AGB stars fall below our observational 
limit. 

In Figure 13 we show the IRAC [3.6] — [4.5] vs. [5.8] — 
]8.0] colour-colour diagram of all the sources in our cata- 
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logue. The extreme AGB stars in our sample are indicated 
by red stars. These sources are often obscured in the optical 
and near-IR due to circumstellar dust extinction {A^ « 3—4) 
and are bright at 8 fim due to thermal emission by the same 
dust. The red colours for the majority of the other sources 
are consistent with mass-losing AGB stars. 

4.4 Mass-loss rates of the AGB candidates 

Using mid-IR colours as a proxy for dust-production rates, 
we can estimate the total dust input from the AGB can- 
didates in M32. We derived mass-loss rates using the 
dependence on [3.6]-[8.0] colour derived empirically by 
Matsuura et al. (2009), for the 160 AGB candidates in our 
catalogue which have an IR-excess of [3.6] — [8.0] > 0.5 
mag. This method of using only mid-IR colours to determine 
mass-loss rates for AGB stars relies on a number of simpli- 
fying assumptions; as such these values should be treated as 
a first-order estimate into the dust production in M32. We 
note that the prescription used to estimate the mass-loss 
rates assumes that the stars in our catalogue have a single 
dust composition, luminosity, wind speed, effective temper- 
ature (Tefi) and solar metallicity. As the colour excess only 
provides a measurement of the dust content, a Galactic gas- 
to-dust ratio of 200 was adopted to convert the measured 
dust mass to a gas mass-loss rate. This conversion factor 
is highly uncertain and it is unclear if this ratio is applica- 
ble to extra-Galactic evolved stars. Due to the large number 
of assumptions we have made about a source, the individual 
mass-loss rates for the candidate AGB stars have high uncer- 
tainties. However, for statistically large samples cumulative 
mass-loss rates should provide a reasonable approximation. 

Individual stellar mass-loss rates range from 6 x 10“® 
to 2.6 X 10“® Mq yr“^ and the cumulative mass-loss rate 
is (1.5 — 1.8) X 10~^ M 0 yr“^, depending on the fraction 
of sources that belong to the evolved stellar population. 
Background galaxies and other interlopers that occupy the 
same region of the GMDs as dusty evolved stars may result 
in some overestimates in the total dust input from sources 
in our catalogue. Dust production is dominated by the five 
most extreme sources in our catalogue, which produce over 
30 per cent of the total dust input. The reddest source in our 
sample contributes 15 per cent of the total. However, this 
source falls outside the colour-cut for x-AGB stars and has 
a strong 24-/rm excess; it is probably an unresolved back- 
ground galaxy, thus we exclude it from our total. 

As the candidate x-AGB stars are centrally concen- 
trated about M32 rather than distributed evenly across the 
IRAC held of view they are less likely to be background 
galaxies. The colour- derived, integrated mass-loss rate for 
the x-AGB candidates is 1.45 x 10“^ Mq yr“^, which corre- 
sponds to 97 per cent of the total dust input in the galaxy. 
The remaining 3 per cent is returned from AGB candidates 
with moderate IR excess. For comparison, we also estimate 
the cumulative mass-loss rate using the colour-MLR rela- 
tion derived by Gullieuszik et al. (2012); this was found to 
be 2.1 X 10“'‘ Moyr-^ 

Due to the severe stellar crowding in the core region of 
M32, not all the red AGB stars with extremely high mass- 
loss rates will be accounted for due to the confusion. Thus 
the cumulative mass-loss for M32 is likely to be a lower 
limit. Other evolved sources in our catalogue produce only 


moderate amounts of dust and gas (less than 10 ®M 0 yr ^) 
and contribute little to the integrated value. 

Overall, our hndings for the total dust input agree with 
other recent studies. In the LMC {Z = O. 5 .Z 0 ) the global 
dust injection rate from evolved stars was found to be on the 
order of 2 x 10~® M 0 yr“^, and for the SMC {Z = O. 25 Z 0 ) 
evolved stars produce 9x 10“^ M 0 yr“^ of dust (Riebel et al. 
2012; Boyer et al. 2012). In both galaxies, extreme AGB 
stars are the primary contributors to the total dust produc- 
tion, accounting for approximately 65 per cent of the dust 
production but only 3 per cent of the sample by number. 
Similarly, the dust input for several more metal-poor dwarf 
irregular galaxies in the Local Group was found to be gov- 
erned by a few AGB stars undergoing intense mass-loss via 
superwind (Boyer et al. 2009). Here integrated total mass- 
loss rates range from 4.4 x 10“® to 1.4 x 10“® M 0 yr“^. 

One notable difference between M32, the Magellanic 
Clouds and the dwarf irregular galaxies is the metallicity of 
the host population. In metal-poor environments, the dust 
yield from oxygen-rich AGB stars is reduced, unlike car- 
bon stars where dust yields do not depend signihcantly on 
the initial metallicity of the progenitor (Groenewegen et al. 
2007; Sloan et al. 2008, 2012). Consequently the fractional 
contribution from oxygen-rich dust may be higher in metal- 
rich environments like M32. 

Javadi et al. (2013) estimate dust-production rates for 
~110 variable stars detected at 8 jj.m in the central square 
kiloparsec of M33 that are of near-solar metallicity. Despite 
there being a comparable number of stars detected at [ 8 . 0 ] 
in M33 as M32, the total mass-return in M33 is much higher 
due to the contribution of several massive, very dusty stars 
with mass-loss rates of a few xl0“® M 0 yr“^. This estimate 
also includes red giant variables for which no mid-IR coun- 
terpart was identified. In contrast to the Magellanic Clouds, 
the dust returned in the central regions of M33 from low- 
mass oxygen-rich AGB stars appears to be comparable to 
the contribution from the more massive carbon stars and red 
supergiants. However, due to the distance of M33 (which is 
not much different to that of M32) and the relatively sparse 
points that comprise the spectral energy distribution, there 
is a great deal of uncertainty in the classification into carbon 
and oxygen rich stars. Furthermore, like M32, the AGB pop- 
ulation detected by Spitzer in M33 is dominated by candi- 
date x-AGB stars with only a small fraction of non-extremes; 
thus there is sufficient information at this time to constrain 
which species of dust is present. 

If a large fraction of the x-AGB stars in M32 turn out to 
be oxygen-rich then the dust return could be several times 
higher; this is due to the lower specihc opacity of silicates 
compared to amorphous carbon grains. Owing to the higher 
metallicity of M32, O-rich stars might be expected to dom- 
inate the dust input at lower luminosities; unfortunately, 
the high surface brightness, limited sensitivity and resolu- 
tion conspires against dust-injection rates being calculated 
for the less luminous (oxygen-rich) sources. Gurrent esti- 
mates for the dust yield from non-extreme sources which 
comprise 25 per cent of stars with an IR excess, show that 
they return approximately 3 per cent of the measurable dust 
production. Hence, AGB stars significantly fainter than our 
detection limits are not expected to show any appreciable 
contribution to the dust-production rate in M32. 
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4.5 Comparison to Davidge (2014) 

During the refereeing process, the publication of an indepen- 
dent reduction of this data-set was made public (Davidge 
2014). Comparisons of the mid-IR CMDs show sources ex- 
tracted by Davidge (2014) reach fainter magnitudes at than 
sources in our catalogue; this is due to differences in the 
level above background chosen for extraction and the strin- 
gent criteria we require for inclusion in our catalogue. This 
invariably means that some real point-sources are rejected 
from our sample, but contamination from blends, non-stellar 
objects, and background galaxies are reduced, enabling a 
more accurate estimate of the dust-production rate for M32. 
Davidge (2014) supplement the mid-IR IRAC data with 
near-IR ground based observations to compare the present- 
day stellar population of M32 with the star-forming history 
measured by Monachesi et al. (2011). Although our meth- 
ods of analysis differ we obtain similar estimates in age (on 
the order of a few gigayears) for the dusty evolved stars in 
M32. 


5 SUMMARY AND CONCLUSIONS 

We have presented mid-infrared observations of the dwarf el- 
liptical galaxy M32, obtained with the Infrared Array Cam- 
era on board the Spitzer Space Telescope. These images re- 
solve individual stars in the bulge of M32, revealing a rich 
population of dusty evolved stars. In Section 3.2 it is shown 
that this population is more centrally concentrated than the 
underlying optical sources, indicating a radial population 
gradient. 

Despite the strong levels of contamination from M31, 
we find that luminous stars with significant dust emission 
dominate the 5.8- and 8-pm luminosity functions. We es- 
timate the tip of the AGB to be M3.6 ~ —10.8 mag and 
resolve stars to depths ~3 mag fainter than this limit. We 
do not reach the red giant branch tip. These luminous AGB 
stars may be the progeny of an intermediate-age population 
in M32 with lifetimes between 0.2-5 Gyr. 

Using mid-IR photometric criteria, we identify 110 ex- 
treme (x-AGB) star candidates, corresponding to approx- 
imately half of all sources detected at both 3.6 and 8 fim. 
These red stars are highly enshrouded by dust, thus they will 
often be missed in optical and near-IR surveys. In the Mag- 
ellanic Clouds the majority of the current dust input into 
the ISM comes from the extreme AGB stars even though 
they comprise < 3 per cent of the population; we expect a 
similar scenario to occur in M32. Using an assumed gas-to- 
dust ratio of 200 we estimate the total mass-loss rate from 
all the x-AGB candidates in M32 to be 1.45 x 10“^ M© yr“^ 
and measure a cumulative dust-injection rate of 7.25 x 10~^ 
M© yr"^ 
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